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Survey 


Four working points are treated in this report. 

1. An inverse integral procedure for calculating separated, turbulent 
boundary layers at infinitely long, crabbing wings is described. 
The concept of calculating inversely into the separation area 
with the boundary layer equations has been extended to flows at 
infinitely long, crabbing wings for the first time here. 

The displacement thickness distribution of the velocity profile 
in the direction of wing depth is utilized as input, and the 
corresponding velocity at the outside edge of the boundary layer 
calculated according to size and direction. It is assumed that 
the velocity profile in the direction of main flow can be 
represented by the expanded Coles profile and the velocity pro- 
file in the direction of wingspan by Coles profile for the case 
of pressure gradient of zero. 

A comparison with measurements in a separating, turbulent 
boundary layer is demonstrated. 

2. An integral procedure developed at Dornier for calculating 
three-dimensional, turbulent, incompressible boundary layers 
has been expanded for the adiabatic case to compressible flows. 

The integral quantities are converted by the Stewartson- 
Illingworth transformation and the calculation carried out in a 
transformed space. As a conclusion of the calculation, the 
calculated quantities are calculated back to the physical space. 
Comparisons with measurements and other procedures are shown. 

3. Calculations results with three-dimensional, transsonic wings 
are shown, taking friction into consideration. The pressure 
distribution is determined iteratively with a potential and 
theoretical (small perturbation) procedure and a three-dimensional 
boundary layer procedure. 

For the case of a turbulent separation of the boundary layer in a 
wing section, an approximation procedure is described for calcu- 
lating the separated, turbulent boundary layer. 

4. The laminar boundary layer is calculated with the Dornier 
integral procedure with tipped rotation ellipsoids and compared 
with results of different procedures and measurements (DFVLR) . 

The three-dimensional displacement surface and the lines of 
turbulent departure are shown. 


^Numbers in the margin indicate pagination in the foreign text. 
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Inverse Integral Procedure for Calculating Separated, Turbulent 
Boundary Layers with an Infinitely Long ^ Crabbing Wing 


Symbols 
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Component of the coefficient of friction in the x 
direction 

Component of the coefficient of friction in the s 
direction 

Component of the coefficient of friction in the y 
direction 

Entrainment coefficient 

Entrainment coefficient in balanced boundary layers 

Form parameter of the velocity profile in the 
s direction 

Form parameter of the velocity profile in the 
X direction 

Form parameter of the velocity profile in the 
y direction 

Form parameter of the velocity profile in the 
s direction 

Form parameter of the outside layer of the velocity 
profile in the s direction for which U>. 0. 

von Kaman constant 

Direction of transverse flow 

Quantities defined in Appendix C 

Direction of main flow 

Component of velocity in the x direction 

Component of velocity in the s direction 

Component of velocity in the y direction 

Component of velocity in the n direction 
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Surface distance to the point in the boundary layer 
for which U = 0 . 
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Angle between the s and x directions 
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Surface flow line angle, angle between the resulting 
surface thrust direction and the s direction 
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Appendix C i e 
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Displacement thickness of the three-dimensional 
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the X direction 
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Momentum loss thickness of the outer layer of the 
velocity profile in the s direction for which U ^ 0 

Momentum loss thickness 


6 



Momentum loss thickness 


Momentum loss thickness of the velocity profile in 
the n direction 


6 

y 


O 

f 

J P. 


o 

6 


e e 


Pv 


P V 
e e 


(1 - ^ dz 
c 



Momentum loss thickness of the velocity profile in 
the X direction 


Momentum loss thickness of the velocity profile in 
the y direction 


6 , e 

a 0 ) u 


Derivations of e according to 6 , n f w and u^ 


V 


Dynamic viscosity 


/6 


n 


Coles parameter 




Surface friction parameter in the y direction 





hr h' h 


u 



1/2 

■) 


Quantities defined in Appendix C 
Quantities defined in Appendix C 
Derivations of according to 5 , w and u^ 

Surface friction parameter in the s direction 


Index 

e At the outside edge of the boundary layer 


4 


1.1 


Introduction 


In spite of the progress in computer technology and in numerical 
methods, numerical solutions of the Navier-Stokes equations for separated 
flows are still very complicated today and not required for a portion 
of the applications. It appears sensible, on the basis of work in the 
past years, to calculate flows, demonstrating separation in a limited 
area, with the boundary layer concept. 

The solutions of the boundary layer equations in the two-dimen- 
sional case, however, become singular when given the velocity distribu- 
tion at the outside edge of the boundary layer (direct problem) in 
the separation point, where the surface friction becomes zero [1,2], 

In the case of the direct problem, it becomes impossible to integrate 
the boundary layer equations downstream from the point of separation. 

Gather all and Mangier [3] have demonstrated numerically, that 
separation cpccurs without leading to singular solutions, when a boundary 
layer quantity, e.g. the displacement thickness or the surface shearing 
stress, is given and the corresponding external flow is calculated 
(inverse problem). Moreover, it is shown in [3], that the usual inte- 
gration direction directed downstream can be retained downstream from 
the point of separation. 

In recent years, several inverse procedures, integral procedures 
[4-9] and difference procedures [10,11,12] for two dimensional, laminar 
and turbulent boundary layers have been developed. In the present work, 
an attempt is made for the first time to apply the concept of the inverse 
procedure to flows with infinitely long, crabbing wings. 

Of course, flows can no longer be calculated with this method, /8 
in which large interactions occur between the flow free of friction 
and that subject to friction, and the boundary layer equations no longer 
have any validity. 


1.2 Fundamentals 


The coordinate directions x and y are shown in Figure 1 with an 
infinitely long, crabbing wing. In addition, the direction of the 
resulting velocity U , supplying the direction of main flow s, and 
the velocity components in the x direction u and in the y direction 
Vg are shown. The angle between the x direction and the resulting 
velocity U is a. The angle between the main flow direction s and 
the direction of the resulting surface shearing force, describing the 
direction of the surface flow line is 3. 


1.2.1 Description of the Velocity Profile Families 

Profile families for the velocity are described for the direction 
of main flow s by the expanded, double-parameter representation of 
Coles [12]. 


5 



( 1 ) 




z l(i)| U 
e ^ 

V 


kLn. 

0) k 


[1 


- cos 


(TT |)] + 5.1 


where k is the von Kamin constant. 

This manner of writing makes possible the representation of 
reverse flow velocity profiles, similar to [4,6] , for separated flows, 

Cf^ < 0 - 


It applies that 
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At the outside edge of the boundary layer 6 it results from equation 1 
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(3) 


Equations 1 and 3 combined provide the representation of the velocity 
profile, employed for calculating the integral quantities 



1 “ ^ [n {1 + cos 
X 0) 


(Trf) -inf] 


(4) 


The profile families in the direction of transverse flow n, 
normally running to the direction of the main flow s, was not formu- 
lated by the empirical approaches of Mager [13] or Johnston [14] , as 
they are employed in the works for calculating three-dimensional, 
turbulent, adjacent boundary flows [15-17]. 

Instead, it is assumed that the velocity profile in the y direc- 
tion, in which the pressure gradient is identical to z^ro, can be 
represented by profiles at a flat plate. With knowledge of the angle 
a and the profile U/Ug in the direction of main flow, the profiles 

in the n direction V/U can be determined. 

e 

Corresponding to equation 4, the following formulation is 
possible : 


^ - 1 - ^ [B {1 + cos (TT f)} - In |] (5) 


6 



For the flat plate, A = 0.0332 and B = 0.625. A and B are /lO 

constant values in each calculation point, supplying a velocity pro- 
file v/Vg, similar to the 1/7 law of exponents. 

It applies that 
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V . 

~ sin o 
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e , e 


( 6 ) 


Therefore it results that 


V ^ sin a r V U , 

U ** cos a * U 

e e e 


or 


u" “ (ir f)Hn |u - A B} - In f {OJ - A>] 
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( 7 ) 


The velocity profiles in the s direction U/Uq, equation 4 and 
in the n direction V/Ug, equation 7, are employed for determining the 
corresponding integral quantities, in the s and n direction. The 
integral quantities in the x and y direction, in which the boundary 
layer development is calculated, can then be determined over the 
angle a in a simple manner. 


1.2.2 Fundamental Equations 

For solving the problem, the following equations are employed: 
X impulse integral equation 


•dx 


+ 0u 
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du 
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dx 



du 
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dx 



( 8 ) 
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In Appendix A, the integral quantities 0^, and A,* are defined /II 
and the calculation of 0^^, and A,* is simultaneously given from 
integral quantities in tne s,n coordinate system. The components of 
the surface shearing force coefficient in the x direction Cfj, is 
explained in Appendix C. 

Equation 3 leads to an implicit expression for w 
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(9) 


Equation 10 differentiated according to s provides 


+ 25 ^ + 6 1 ] ±Ji; 

ds w ds ‘‘or l(u] ds 


0 e 

U ds 
- e 


( 10 ) 


wing) 


For a random quantity Q, dQ/dy e 0 (infinitely long, crabbing 


ds • U dx ' 
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( 11 ) 


The surface friction equation results from this for equation 9 


~ + 25 7-^+6 [V + A-] 

dx « dx ‘‘or joj[ dx 


. dU 
0 e 

U dx 
e 


( 12 ) 


Since entrainment occurs in the outer area of the boundary layer 
and the transverse flow velocities are small there, the two-dimensional 
manner of consideration can be extended to the flow at infinitely long, 
crabbing wings, similar to the integral procedure [15-171. The 
entrainment is coupled in this case to the velocity profile in the 
direction of main flow. 
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The entrainment equation is expressed as: 


/12 
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= F 


( 13 ) 


The lag-entrainment concept of Horton [l8l is employed for calcu- 
lating the entrainment coefficient F in equation 13, taking into 
consideration the so-called upstream history effects on the coefficient 
F in a non-balanced boundary layer. 
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( 16 ) 


In the case of adjacent flows 
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On 


011 und Hi =t Hi. 


In the case of a separated flow in the direction of main flow, 
the entrainment is related to integral quantities in a similar manner 
to [4], formed from the velocity profile situated above the U = 0 line. 

For separated f lo ws , t he calcul.ation of the U = 0 line and the corres- 
ponding values <Si and <Si * and dn is given in Appendix B. 

With equation 11, it follows for equation 13 /13 
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An approach similar to that for in equation 10 is utilized 

for calculating the surface shearing fofce component Cfy. 

With Cfy = 2 it follows that 


do 

dx 


-0.03243 i ^ 
0 dx 


(18) 


Equation 18 states that the deviation in a, i.e. c„ , is coupled 
to the variation of 6. A prerequisite in this case is tliat the 
boundary layer thicknesses of the velocity layer profiles are identi- ( 

cal in the s, n, x and y directions. 


1.2.3 Solution of the Resulting Equation System 

Five equations are available for the dependent variables n, to, 
Ug, a and P 

X momentum equation 
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=’n ST 
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dx 
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dx 
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entrainment equation 
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lag- entrainment equation 
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surface friction equation in the s direction 


. dH 


do) 
0) dx 


+ P3 


du 
e 

dx 


- 6 


6 * 
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d 6* 

> 

dx 


( 24 ) 


surface friction equation in the y direction 


da 

dx 


- 0,03243 


1 ^ 

6 dx 


(25) 


The coefficients resulting in equations 21, 22 and 24 and the 
expression d6/dx are explained in Appendix C. 

The equations 21 - 25 are solved with an explicit intermediate 
step procedure. As initial values, values for 0xi, H , 3, and u are 
inserted and the distribution of the displacement thickness is 
prescribed as a function of x. 


1. 3 Comparison with Measurements 

Van den Berg and Elsenaar fl9] have measured the development 
of a turbulent boundary layer on a crabbing plate (sweep angle 35°) 
in the range of low velocity. The wind tunnel wall above the plate 
was designed in such a manner that the generated pressure increase 
was large enough to cause boundary layer separation. Moreover, the 
attempt was made to approximate as closely as possible the conditions 
of the infinitely long, crabbing wing by the form of the end plates 
on the measurement plate. 

The measured velocity profiles v/v^ are plotted in Figure 2 /15 

over the surface distance z standardized with the boundary layer 
thickness 6. Flow separation was observed near the measuring point 
8 in the experiment. As can be seen in Figure 2, the requirement 
made for the calculation procedure that the velocity profile in the 
y direction can be approximated by the velocity profiles at the 
flat plate (drawn line) is confirmed here by the measurement results. 

In the adjacent, separating and separated range, the v/vq profiles 
are sufficiently well represented by profiles at the flat plate. 

The upper diagram in Figure 3 shows the measured quantities of 
6 *, employed as input for the inverse procedure (the drawn line is 
tne curve fit employed for the input) . In the lower portion, the 
displacement thickness of the main flow profile 6,* has been compared 
with the measurements. The calculation results are a good representa- 
tion of the measured quantities 6^^*. A comparable good quality of 
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calculation for the quantities 6 iif 0 i 2 r 622 62 * is given in the 

Figures 4 and 5. Figure 6 shows the boundary layer thickness 6 
and the surface flow line angle 3 . 

The form parameter H (relationship of displacement thickness to 
the momentum loss thickness) of the velocity profiles in the direction 
of main flow Hg, in the x direction Hjj and in the y direction Hy are 
shown in Figure 7. The assumption is again confirmed here that there 
is a plate boundary layer velocity profile in the y direction. For 
this case, Hy = 1.29. For the entire area, the measurement also shows 
a constant value for Hy, corresponding approximately to the said value 
of 1.29. The comparison of measurement and calculation of the values 
for Hg and H^ is very satisfactory. 

The measured and calculated values of the surface shearing force /16 
components c^j^, Cfy and Cfg are shown. The value cfx = 0 indicates 
separation at an infinitely long, crabbing wing by definition. As can 
be seen, the calculated separation point is situated in the area of the 
experimentally observed separation. 

In Figure 9, the finally desired results of the inverse boundary 
layer calculation for the external velocity component u^ and the angle 
a have been compared with the measured quantities. Three different 
interpretation possibilities for the measured external velocity data 
can be given: 1 . Ug and a measured, 2 . u^ measured and a calculated 

from the condition of infinitely long, crabbing wing and 3 . Ug and a 
calculated from the measured surface pressure distribution with the 
condition of infinitely long, crabbing wing. As can be seen, slight 
differences between the different representation possibilities only 
result near the separation point. The calculation provides a good 
agreement for the values of u and o. The calculation results are 
very closely situated to the quantities, derived from the measured wall 
thickness . 


2 . Expansion of the Integral Method on the Calculation of /17 

Three-Dimensional, Turbulent Boundary Layers to Adiabatic , 
Compressible Flows 

Symbols 

a Speed of sound 

Cj Transformed surface coefficient of friction 

F Transformed entrainment coefficient 

M Mach number 

T Static temperature 

T* Reference temperature 

z Surface distance in the physical space 
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z 

Y 
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Surface in the transformed space 
Ratio of the specific heats 
Transformed boundary layer thickness 





Transformed displacement thickness of the velocity 
profiles in the direction of main flow 


5 

5* = - [ — dz Transformed displacement thickness of the velocity 

i profile in the direction of transverse flow 


\ Transformed momentum loss thickness of the velocity f 

01 j ts I (1 _ IL) dz profile in the direction of main flow 

Index /18 

e At the outer edge of the boundary layer 


o 


Reference state 


2.1 Introduction 


In the reference [ 20] there is a calculation method developed at 
Dornier described for three-dimensional, incompressible, turbulent 
boundary layers. This procedure is a further development and improve- 
ment of an existing procedure developed by My ring 1 2l] . The differences 
to the work of Myring are: 

1. The single parameter exponent profile for describing the velocity 
distribution in the direction of main flow has been replaced by 
double-parameter Cole profiles [12 1. The surface friction para- 
meter describes the velocity distribution near the surface and 
the pressure gradient parameter determines the velocity distri- 
bution in the external section of the boundary layer. 

2. The empirical information on the surface coefficient of friction 
additionally required in using exponential profiles is not 
required . 

3. Instead of the entrainment method employed by Myring in which 
the entrainment coefficient is determined directly as a function 
of the form parameter, a "lag entrainment" method has been 
utilized in the present method, taking into consideration the 
unbalance effects of the boundary layer. 

In the following, the alterations are described for expanding the /19 
procedure to the calculation of adiabatic, compressible flows. 
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2.2 


Fundamentals 


The effects of compressibility on the boundary layer development 
are taken into consideration by the Stewartson-Illingwor th transfor- 
mation 


dZ « 


£_ 



( 26 ) 


The calculation procedure is limited to adiabatic flows and mach numbers 
M ^ 2. With the example of the transformed momentum loss thickness of 
the main flow profile 


'll 


u ,, u 
^ (1 - dz 
e o 


in which the thickness p no longer appears, the relationship between 
the compressible and the transformed, incompressible quantities is 
demonstrated . 

Utilizing equation 26, it can be shown that 
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( 27 ) 
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2 ^ Y-l ^ 
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and ~ = (1 

2 e , . 


Equation 27 can also be written in the same manner for the momen- 
tum loss thicknesses 


5i2 •» C 012 

S 21 *> C 02 1 (28) 

622 = 5 02 2 , 


and for the displacement thickness of the transverse flow profile 


_» * 
62 “ C ^2 


When the quadratic temperature-velocity ratio and the condition 
dp/dz = 0 is utilized, the following formulation can be made: 

^ i + 2 [1 . (^)2 _ (^) 2 j ( 30 ) 


With equations 26 and 30, the result for the displacement thick- /21 
ness of the main flow profile 


61 = C * l 5 l + (Si +611+ §22) 1 , 
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and for the boundary layer thickness 


6 = C * (3 1 + §11+ §22)3 


( 32 ) 


The coefficient of friction cj^^ can be written with the aid of 
the Eckert reference temperature concept 


- Te 
®fx “ yX 

with 
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— =■ 1 + 0.13 M * 

-.X e 


(33) 


(34) 
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The entrainment coefficient F can be given in the same manner 
according to a proposal by Horton [18] 


F 



(35) 


In the derivation of the boundary layer equations in the trans- 
formed space, the following conversion has been utilized 
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with similar expressions for 
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With the exception that the velocity profiles in the transformed /22 
states can be represented by those of the incompressible flow, the 
resulting equations are expressed as: 
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X- impulse equation 
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surface friction equation 
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entrainment equation: 
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lag-entrainment equation; 
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with 
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The three-dimensional displacement thickness S * in the physical /24 
space is calculated with 
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The quantities occurring in equations 36-42 are defined in 
reference [20 1. The desired quantities in the physical space, e.g. 
011 , are determined with the aid of equation 3 from the calculated 
quantities in the transformed space, 9 ii. 
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2.3 


Results 


The boundary layer measurements in a supersonic jet carried out 
by Hall and Dickens 1221 have been recalculated. The shape of the 
supersonic jet, the course of the flow lines A, B and C and the Mach 
number distribution along the flow line B are presented in Figure 10. 

Figure 11 presents the results of the quantities 6 ii, H and 3 . A 
comparison with the calculated results of the procedure by P. D. Smith 
l23J is also indicated. 

A further calculation result of the procedure described here 
has already been presented in the report of the workshop on three- 
dimensional turbulent boundary layers [24] . Figure 12 shows the 
foundation shape of the wing. Figures 13 and 14 present the pressure 
distributions for a flow approach mach number of M =0.5 and an angle 
of inclination of 0 or 8 *^. Figures 15 and 16 present the calculated 
results of the present procedure (complete circles) compared to the ( 

average result of all calculation procedures participating in the work- 
shop (three integral procedures and five difference procedures) . 


3. Pressure Distribution Calculations with Three-Dimensional , /25 

Transsonic Wings, Taking i nto Consideration the Friction 
Effects 


3. 1 Introduction 

The alteration work in boundary layer procedures for the approxi- 
mate calculation of separation areas have been undertaken within the 
framework of the present commission. Results of this work have also 
been presented in the ZKP report, Dornier Report No. 79/38B. 

In order to calculate the frictionless flow around three-dimen- 
sional wings in the transsonic range, the small perturbation procedure 
[ 25 ] is utilized. The development of the three-dimensional boundary 
layer is determined with an integral procedure [20,26] . Iterative 
calculations of the frictionless flow and flow subjected to friction, 
using the "displacement surface concept", produce the pressure distri- 
bution in the converged state while taking into consideration the 
effects of friction. 


3 . 2 Approximate Determination of Boundary Layer Quantities in 

the Case of Separation 

When laminar separation occurs in the calculation before achieving 
the defined conversion point, the calculation is continued with the 
condition of an adjacent, turbulent boundary layer. In the case of the 
turbulent boundary layer separation, no further data can be determined 
downstream from the separation point with the available method. 

To avoid interruption of the cycle of the iterative, frictionless /26 
and friction calculation, the following, approximate determination of 
is applied, see Figure 17. 
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If a turbulent separation occurs in the calculation in the wing 
section 7^+4 ^t the percent line the wing section Yn+4 is excluded 

from the normal calculation. In the wing sections y^+s + Yn+S' 
boundary conditions for the remaining areas of the three-dimensional 
calculation are determined by assuming the condition of the infinitely 
long, crabbing wing in these wing sections. Finally, after further 
integration in the x direction, separation in the wing sections yr. 4.9 t 
Y n +6 determined. " ^ 

These wing sections are now excluded from the normal calculation 
and boundary conditions for the remaining areas of the three-dimensional 
calculation are determined in the wing sections Yn+i snd y„ -y. In 
this manner, a separation line can be calculated in the case of three- 
dimensional wing boundary layers. In the wing sections excluded from 
the normal calculation, the boundary layer development is determined 
in an approximation. In order to calculate the desired displacement 
thickness in this case, only the entrainment equation is employed with 
a defined form parameter of the boundary layer, corresponding to a 
separated turbulent boundary layer. Moreover, the value of the surface 
flow line angle determined at the separation point is maintained. For 
this reason, the solution of the x and y-impulse integral equations 
is not required. 


3.3 Example Calculations with Wings /27 


3.3.1 PT7 Wing 

The fundamental shape of the PT7 wing, developed and measured by 
the FFA in Stockholm, is shown in Figure 18. The results of the cal- 
culations are compared in wing sections 1 -f 3 with measurements at a 
mach number of 0.9 and a Reynolds number of 1.38 • lO^ [l/m]. 

Figure 19 shows the three-dimensional displacement thickness in wing 
section 2 related to the root wing depth, plotted over the relative 
wing depth x/C for a pure turbulent flow. Curve 1 corresponds to the 
boundary layer calculation with the first frictionless pressure dis- 
tribution. The results are smoothed and only 15% of the determined 
displacement thickness are transmitted to the second frictionless 
calculation. The results of the second boundary layer calculation, 
curve 2 , are again smoothed and 80% of the displacement thickness 
transmitted. After the third boundary layer calculation, smoothing is 
no longer carried out and 100 % of the displacement thickness are 
employed. The fifth subsequent calculation did not demonstrate any 
deviations compared to the fourth, and the result was considered as 
converged. 

As can be seen, turbulent separation occurs in the first boundary 
layer calculation at a 85% wing depth. The values of 6 * for x/c > .85 
correspond to the above-described calculated approximation. In the 
second boundary layer calculation, the frictionless pressure distribu- 
tion was altered in such a manner that separation no longer occurred. 
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For, comparison, the converged result for a boundar 
IS drawn (dotted line), for which up to 20% wing depth 
flow was assumed. 


y 

a 


layer flow 
laminar 


‘he pressure distributions determined 
by calculation in the wing sections 1 ^ 3 in comparison to the 
measurements. In all wing sections, a marked upstream motion of the 
area of pressure increase is determined on the upper side of the winq 

‘he pressure dlstriLtion remains ataSs; 
front of the area of pressure increase, only in the wing 
mined^^ ^ ^ clear reduction in the under pressure level be deter- 


The rear edge pressure on the lower side of the wing is reduced 
in all sections. 

. . . ^^ 9 ure 23 shows the comparison of the converged pressure distri- ^ 

bution in the case of pure turbulence and a laminar-turbulent boundary 

turbulent case, the area of pressure increase 
IS shifted more greatly upstream, explained by the fact that the inter- • 
between boundary layer and external flow becomes more intensive 
in the case of larger boundary layer thicknesses. 

Figure 24 shows the separation areas on the wing, resulting in the 
course of the iterations. It becomes clear that the separation area 
has become very much smaller in the converged case than in the first 
subsequent calculation. This result underlines the practicality of 
the approximate determination of boundary layer in the "separated" area. 


3.3.2 ZKP Winq 

Figure 25 shows the basic form of the wing for the ZKP winq. 

The displacement thickness for the wing section 2 is shown in 
Figure 26 in the course of the iterative calculations. Figures 27-29 
show the measured pressure distributions and those determined in 
calculations. The comparison shows clear deviations between measure- 
ment and calculation, partially attributable to elastic shaping of 

the wing and perhaps to the course grid employed in the frictionless 
case. 
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Laminar Boundary Layers in Tipped Rotation Ellipsoids 


4 . 1 Introduction 

A detailed description of this work can be found in reference [27], 
Tne calculation results shown here are determined with the Dornier 
integral procedure for calculating three-dimensional, laminar boundary 
layers [26]. ^ 
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4.2 


Results 


The frictionless external flow is determined with a potential, 
theoretical method. The lateral ratio of the observed ellipsoid was 
6:1 and the angle of inclination a = 10*^. 

The calculation results are compared in Figure 30 with those of 
a difference procedure [28]. 0ii and 6i* are the momentum loss thick- 
nesses or displacement thickness of the velocity profile in the direc- 
tion of the main flow, T is the resulting surface shearing force and 
3 the surface streamline angle. 

As Figure 30 shows, the results of the integral procedure are in 
good agreement with those of the difference procedure. 

The Figures 31 and 32 show results for the surface shearing force 
\ in comparison with measurements [29] and results of difference 
procedures . 

The development of the three-dimensional displacement thickness /31 
6* is shown in the Figures 33 and 34. Figure 35 describes the course 
of the surface streamlines. The vectors drawn in indicate the direc- 
tion of the resulting surface shearing force. The calculated turbu- 
lence departure lines are shown in the Figures 36—38 arid compared 
with the results of the difference procedure. 


5. Summary /32 


5 . 1 Inverse Integral Procedure 

An inverse boundary layer integral procedure has been developed 
for calculating separated, turbulent boundary layers with infinitely 
long, crabbing wings. The displacement thickness distribution of the 
velocity profile in the direction perpendicular to the leading edge 
of the wing is employed as input. The velocity distribution is cal- 
culated according to size and direction at the outer edge of the boundary 
layer . 

The expanded double-parameter profile family of Coles is employed 
for describing the velocity profile in the direction of wing flow. 

It is assumed that the velocity profile v/vg in a direction parallel 
to the leading edge of the wing corresponds to that on a flat plate. 

The entrainment equation is employed with a lag-entrainment 
method for calculating the entrainment coefficient, the surface fric- 
tion equation is employed for determining the surface shearing force 
component in the direction of main flow and a simple relationship is 
employed for determining the surface shearing force component in the 
direction parallel to the leading edge of the wing in order to solve 
the problem in the x-impulse integral equation. 

A comparison with measurements in a separated, turbulent boundary 
layer confirms the assumption made on the v/Vq velocity profile. 
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Moreover, the quality of agreement demonstrates the usefulness of the 
calculation procedure. 


5*2 Expansion of the Three-Dimensional, Turbulent Integral /33 

Procedure to Compressible Flows 

The expansion of the existing integral procedure for three- 
dimensional, turbulent flows to compressible, adiabatic flows has 
been made with a compressibility transformation. The calculation 
examples show good agreement with measurements and other procedures. 


5.3 Example Calculations with Wings 

The iterative calculation of pressure distribution with three- 
dimensional transsonic wings while taking into consideration the 
effects of friction produce a satisfactory agreement between measure- 
ment and calculation. The usefulness of an approximation procedure 
for determining the boundary layer development in the case of separated, 
turbulent boundary layers has been demonstrated. 


5 . 4 Laminar Boundary Layers with Tipped Ellipsoids 

A comparison of the calculation results of laminar boundary 
layers in the case of tipped rotation ellipsoids of the Dornier 
integral procedure with available measurements (DFVLR) and 
difference procedures produced good agreement. Even the position of 
the turbulent departure lines at the ellipsoid has been calculated 
in very good agreement with the difference procedure. 
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Appendix A 


Calculation of the integral quantities 0n and Ai* from integral 
quantities in the s, n coordinate system. 

It applies that 


u u V 

iT “ ip cos a - — sin a 
e e e 


(la) 


V U ^ V 

u~ “ u” tj~ ® 

e e 


(2a) 


It results that 


©11 =*011 cos^ a - (012 + p2i) sin o cos a + 622 sin* a 

A » - * . « 

Ai => 61 cos 0-62 sin a 


(3a) 

(>ia) 


The expression with equation 4 and 7 is 

6,j => 6 g {[n M + I] - g [1.5 n* + 3.1794 M n + 2 ]} 

012 = tg a ^ {1.5 n* u)* - 0.031163 II Iwj + 3.1794 [II o) |u)l 
- 0.01662 (II |u)l + 0.625 w) ] + 2 (1) (d) - 0.03324)} 

022 = - tg* a ^ {1.5 n* w* + 3.1794 II |o)| o) - 0.16801 
H |b)| + 20)* - 0.19901 U) + 0.00505} 

= 5 ^ tn 1] 

$2 “ “ ^ ® {w[n + 1] - 0.05402} 


(5a) 


(6a) 


(7a) 


(8a) 


(9a) 



For 021 it applies that 


ZM 


02 1 = 012 +62* 


(lOa) 


Appendix B 

Calculation of the entrainment is related to integral quantities, 
similar to the procedure in [ 4 1, formed from the velocity distribution 
situated above the U = 0 line. 

According to equation 4, it applies for the position z^, at 
which U = 0 ^ 


1 o g [n M {1 + cos dr - In 

or 

Z = 6 exp [n {1 + cos <rr f) - ii] 
o 0) 6 0) 


(16) 


(26) 


It follows for the integral quantity with to < 0 



Si 


- 6 




f (1 - ^) (l.'s - 4H + 2 ) - 


3- In T- (H + 1) 


- 3- (In - - sin (TT jS) - sin {2rr 5^) + 2 | (1. 


289 


(TT^) 


3 


In ^ sin dr + 


hr t) 


2 , s 


35280 3265920-’ 




18 


.600 


(5b) 
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Appendix C 


For the components of the surface shearing force coefficient it 
applies that 

Cfx “ ®fs “ si*' ° tan $) (Ic) 

*^fy “ ®fs a + cos a tan 0) (2c) 


Eliminating the value tan 3 from equation Ic and 2c, 


'fs 


'fx 


cos a 


[1 


*^fs 


sin a] 


(3c) 


The following derivations are to be formed for determining the 
coefficients in the equations 21, 22 and 24, where it applies for 
example that aoj,/36 = 


Gii, 


in. 

6 


0I2i 


On 

6 


(4c) 


62 X I 




Ga 1 
6 

« 

in 

6 


6 


_ 6a 

^25 - T 


Gixjj = I |<ol - 6 ^ [3n + 3.1794 


Giarr = tg a ^ {311-w^ + 3.1794 w lojj - o.o84o Iu|} 


Gain “ 

Gaan- ” - tg* ct ^ {3n + 3.1794 |u| (U - 0.16801 la)|> 

11 K 


(5c) 







011 


fw 


u> 


6i2 


(l> 


*tg o ^ l«j| 

(tjl M. + 1] - 2 ^ [1.5 + 3.1794 -y> + 2]} 

k w ^ 

tg o fr {3 n* w - 0.031163 H + 3.1794 [2 lw| - 0.01662 H -y-] 
+ 4 (i)> 


/44 


021 


(I) 


022 


b) 


01 


(I) 


62*' 

*W 


" ®l*0, " 

S -tg^ a ^ {3 CO + 6.3588 H lw| - 0.1681 H + 4 co - 0.19901} 
k 

6 _ ( 6 c) . 

*■ k"V- 

» -tg a jj- n 


For example, with e», = 30,,/3u^ and3d/3U^ = 

It results that u e e 


- sin a cos a I/O 


0iiu ” o 


0 X 2 u 


_ 9x2- 

U 

e 


02 lu 
022 u 


®12u ^2 u 


fiZi. 


61 . 


62' 


o 


62 
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When A is representative for 6 , II and co, the result for the 
coefficients in equations 21 , 22 and 24 is 
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011 . « ^ ” { 012 . ■.+ 02 1.} Sin a cos a + 622 sin* a 


(8c) 




( 9 c) 


For 0 x 1 and A*i the result is 
u u 


01 




11^ cos* a - {0i2^ + 021^} sin a cos a - 622 sin* o 


+ sin* a cos* a — [2 0n + « - sin* g . a \ i o 1 

Ug “ 'sin a cos a ^^“*2 + 821) - 2 $22! 


(loc) 


^‘.u “ cos a - 62^ sin o + sin o cos a ~ [6^ sin a + 62* cos a] 


1 r,* 


(11c) 


The result for the coefficients i|) 
u. 


Ax 

e 


"6 “ ir-Ar^ 

e 


ij; » ^ (2 

ri 


a 


r) - 


e 




n 


n 


% 


0) 


‘to 


(12c) 


and moreover 


2 6 -!iiL 

to 




( 13 c) 
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The quantities are defined as follows: 




0. u 


P2 




Pi 


U 

5 

0 2 
o 


For the derivation 6^^* it applies with 


f # _ 96 


and equations 4a, 8a and 9a 


« 

Ax 

6 

ts •>« 

1 

k 

COS a 

A 

w 0) 

(hM 



0) 


» il 

A 

cos a 


(14c) 


For the derivation of the boundary layer thickness, it can be /47 
written as follows: ^ 


Si 

<lx 


fiAx air 


s — 

0 ) dx 


» dJI 

+ 

V dx 



(15c) 
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with 


cos a ^ {ji Jwl ^ . 

fi / 

A • 

® -k Aj cos o ^ fw| 

# 

5 » ^ cos ct 

“ ^ «e “ 0-10804 cos* a) 


( 16 c) 
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Definition of the Velocity Components and 
the Angle ^ and^ 

a. Direction of the resulting stirface 
shearing force 


•^1 l« 



0 0.2 0.4 0.6 0.8 1.0 1.2 

V/Vg 


Pig, 2: Oomparison of the measured turbulent 

velocity profile v/v^ with the velocity 

profile on a flat plate . [19] 

Key: a. Measurement points no, 
b.flat plate 
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Messung 

geglQttGte Eingobe 


o Messung 

— vorliegendes Yerfahren 


Pig, 3 : Distribution of the measured quantity c * [19] as 


input and the comparison of the calculated quantity 
wdfch the meas\irement [ 19 ] • 


Key; 


a. measurement 
d, calculated 


b, curve fit 
e, solution 


c. present procedure 
f. measured 






lil 

A 


o Messungen 
6 

— vorliegendes Yerfahren 


1 ‘ a 

Y^bercchncle 

b Ablbsung 
Vgemesserte 


X Im] 


Pig, 4: Comparison of the calculated quantities 0n and 9l2 

with the measurements [19] . 

Key; a, calculated h, separation c. measured 

d, measurements e. present procedure 








X Imr 

Pig, 6; Comparison of the calculated boundary layer thickness ^ 
and the sxirface streamline angle with the measurement 

CI9] , 

Key: a. surface d. calculated 

b. measurements e. separation 

c. present procedure measured 





Fig. 7: Comparison of the calculated form parameters 

H , H and H with the measurements [19] . 
s X y 

See Fig. 4 for key. 















Fig, Si Comparison of calculated values of the velocity 

component u^ and of the angle cc with the measured 

values. [19] 

Key: a. calculated d, from the condition of infinitely 

h, separation long, crabbing wing 

c. measured e. surface pressure measurement 










10 15 20 25 30 35 40 45 

X{in) Stromlinie B 


Pig, 11; Comparison of the calculated values of B , IT and 
with measurements [22], 

Key: a, measurement 

h, present procedure 


c , streamline 
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Pig. 13: 


Pressure distribution at the angle of 
inclination oC = o° 





/61 
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Pig, 17; Diagram for explaining the manner of calculation 
in the case of separated flow. 
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e displacement thickness on 
the PT7 wing in the cotirse 






Pig, 20; Pressure distribution at the PT7 wing (section 1) 
in the course of the iterations and comparison 
with measurements. 
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Pig. 21: Pressure distribution at the PT7 wing (section 2) 

in the course of the iterations and comparison 
with measurements; 
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Fig, 32: Comparison of the calculated sircface shearing 

force 2;^ on the line of symmetry facing the wind 

with measurements[29] • 
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Fig, 33: The stirface displacement distribution of an 

ellipsoid with a side ratio of 6:1. 
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Fig. 35: Illustration of the direction of the resulting 

surface shearing force. 
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VEW FROM BELd'.V 



Fig, 36: Comparison of the turhulence departure lines for 

the ellipsoid with a side ratio of 6:1, 






















Fig, 38; Comparison of the turbulence departure lines for 
the ellipsoid with the side ratio of 2:1, 
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